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ABSTRACT 



In this article we prove several important properties of 2*^ p minimum aber- 
ration (MA) designs with k > n/2, where n = is the number of runs. 
We develop a simple method to build MA designs of resolution III. Further- 
more, we introduce a simple relationship, based on product of polynomials, for 
computing their word-length patterns. 
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1 Introduction 

2^~P designs are commonly used in experimentation, where k is the number 
of factors and n = 2''-P is the number of runs or observations. These 
designs are economical, statistically efficient and simple to analyze (see Box 
and Hunter, 1961). 

In many applications, especially in industrial experiments, it is often 
necessary to determine which factors among a large number of candidates 
could affect a particular response. For a given number of factors k, a usual 
procedure is to select a first fractional 2^~p design to accommodate the fac- 
tors with the smallest number of runs. Confusion of effects may lead to 
more than one plausible explanation of data. The initial design should be 
chosen in such a way that simplifies as much as possible the analysis stage 
and reduces the size of subsequent experimentation; see Box, Hunter and 
Hunter (1978) or Tiao (2000) for examples. When the experimenter has lit- 
tle knowledge about the relative sizes of the factorial effects, the minimum 
aberration (MA) criterion selects designs with good overall properties. A 
detailed discussion on the MA criterion can be found in Fries and Hunter 



(1980) and Chen, Sun and Wu (1993). Since Fries and Hunter (1980) intro- 
duced this criterion, many articles have been devoted to find MA designs and 
to study their characterization and structure, for example. Franklin (1984), 
Chen and Wu (1991), Chen (1992) and Tang and Wu (1996). 



In addition to their direct practical application, minimum aberration de- 
signs are also useful in the construction of more complex designs. For exam- 
ple, Ankenman (1999) includes four-level factors in a design, and Bingham 
and Sitter (1999) and Huang, Chen and Voelkel (1998) use them to obtain 
minimum aberration two-level split-plot designs. 

MA designs can be obtained with reduced computational effort, particu- 
larly when the number of runs is low (16, 32 or 64). In this paper we propose 
an original procedure to obtain MA designs of resolution III in a simple way. 
This result could be useful from a practical point of view when the number 
of runs is high. That is not, however, our main interest. Following the work 
of Tang and Wu (1996), we wish to understand the structure of MA de- 
signs, establishing some of their characteristic properties, and to place them 
within the global structure of designs with a fixed number of runs. We will 
see graphically that this structure presents interesting properties. 

In this paper we shall only consider screening designs. Thus, we will 
consider '^-v designs which allow us to study k factors with a minimum 
number of runs n = 2*^"^*, i.e., n/2 < k < n — 1. There are two reasons 
for this choice: first, they arc preferred when the number of factors to be 
analyzed, k, is large; this is precisely the case when the search for a good 
design is more complex. Second, using a symmetry property (Tang and Wu 
, 1996) the properties of designs with k < n/2 can be deduced from those 
with k > n/2. As a consequence, a complete characterization of designs 
having k > n/2 provides a characterization for all designs with n runs. 
Moreover, all designs with k > n/2, except for one, are of resolution /// 
(for k = n/2), and their comparison would require precise criteria to be able 
to discriminate among them. 

Finally, and in addition to providing a characterization of MA designs 
for > n/2, we also provide a way to replace the tedious computation of the 
word-length pattern vector by the computation of an elementary product of 
polynomials. 

This paper is organized in four sections. Section 2 introduces the problem 
and the notation. In Section 3 the main result for this paper is presented, 
namely the necessary condition for any 2'^~p design with A; > ra/2 to be of 
minimum aberration. Sections 4 and 5 are concerned with the computation 
of the word-length pattern, and from the results in these sections we complete 
the conditions which characterize MA designs. 
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2 Definitions and motivation 



A complete two-level factorial design with m factors Bi,B2, contains 
the n = 2"* possible combinations of factor levels. Usually, - and + signs are 
employed to indicate the two levels of each factor and the design is presented 
in a design matrix, which has m columns, one per factor, and 2™ rows corre- 
sponding to the factor treatments. The m columns are independent in that 
none of them can be obtained as a product of the remaining ones. To sim- 
plify the presentation, we shall call these columns or factors the basic ones. 
By multiplying these columns in all possible ways we can generate up to 
s = 2*" — m — 1 additional ones, which will be identified by the combination 
(juxtaposition) of the basic factors which generate them Tri{Bi,B2, 
or more concisely as TTi. The individual elements of column TTj are ob- 
tained by multiplying the corresponding elements of the columns specified 
in Tri{Bi,B2, ...,Bjn)- The collection of the 2"^ — 1 columns 

Hm = {-61,-62, ...,Bm,Tri,TT2, ...,7rs} 

is the key set through which any 2''~p designs with 2"* runs can be defined. 
In coding theory the set Hm is known as a Hamming code. For example, if 
we use A,B,C,D for the four initial columns of the complete 2^ design, we 
obtain 

H4 = {A, B, C, D, AB, AC, AD, BC, BD, CD, ABC, ABD, ACD, BCD, ABCD}. 

Multiplying any column of H^n by itself we obtain the / column, where all 

signs are +. The set {I}UHm with the defined operation has a structure of an 
abclian group, where / is the null clement. The result of the product of sign 
columns can also be derived by looking at the products of the combinations 
of letters which represent them. For example, the product of columns ABC 
and ABD produces the column CD, since AA = I, BB = I. 

Definition 1 A subset d of k columns of Hjn represents a 2*^""^ design. 

For example, the following subset of 9 columns of i?4 

di = {A, B, C, D, ABC, ABD, ACD, BCD, ABCD} (1) 

is a 2^~^ design. 

The advantage of this definition derives from its generality; it enables us 

to define any design with n = 2"^ runs, whether it is a complete factorial 
design or a fractional design 2^~p [k > m). Moreover, through this definition 
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we can make use the tools of set theory. We will then say that a given design 
contains another, is the union of two others, and so forth. 

Many subsets formed by different elements of Hm provide equivalent 
designs which for our work will be considered the same. Two designs are 
equivalent if it is possible to define an isomorphism between them (see for 
example, Tang and Wu, 1996). Thus, Pu (1989) proved that there are only 
five different 2^^^ designs, so that any subset of H/^ with nine elements is 
equivalent or isomorphic to one of those designs. If we consider all designs, 
with any number of factors, the number of different designs with n = 16 
runs is reduced to 45 (46 if we include 0). In figure 1 we present a graph 
showing the relationship between them. Each point or node in the graph is 
a design. The starting point is the empty set (A; = 0); for A; = 1 all designs 
are equivalent and are represented by a single node, and the same occurs for 
k = 2. For k = 3 there are two different designs, {A, B, C} and {A, B, AB}; 
for k = 4 there are three, and so forth up to /c = 15 which is the H4 design. 
The lines connect pairs of designs which only differ in one element, which 
in some cases (solid lines) is indicated on the arc. The design with k + 1 
columns is obtained by adding the column which appears on the arc that 
connects it to its predecessor. The designs for k > 8 can be obtained by 
symmetry, starting from H4 and eliminating the columns specified in the 
symmetric arcs. 

One of the most important properties of this graph is its symmetry, 
due to the fact that if two designs di and d2 are equivalent, so are their 
complementaries di and d2- The notation d will be used to represent the 
complementary of a given design d, where d = Hjn\d. Thus, within the 
subsets of i?4, there are four different designs for A; = 5, there exist four 
others for k = 10, which correspond to the complementaries of the first 
ones. In figure 1 the complementary of a design with k columns is the 
symmetric with 15 — columns. The relationships between complementary 
designs are analyzed by Tang and Wu (1996). In accordance with this, the 
properties of one half of these designs [k < n/2) can be deduced from the 
other half k > n/2. In this paper we shall concentrate in the analysis of MA 
designs for A; > n/2 

The goodness of a design is a consequence of the relations of dependence 
between its columns. Let us consider the design d C H„i, with columns 
d = Ui,?2>->^ik}- A word of length i consists of i elements Cji, ^jj' "•' ^ji 
from d such that 

^n^n--<J.=I (2) 
where I denotes the column having all +'s. The set of all distinct words 
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formed by products involving elements of d gives the defining relation of the 
design. The words in the defining relation correspond to all the interactions 
of the k factors that are confounded with the mean (represented by column 
/). The vector 

W{d) = [ai{d),a2{d),...,ak{d)] 

is called the word-length pattern (WLP) of the design d, where ai{d) is the 
number of words of length i. We shall use in the future instead of ai{d) 
whenever this does not lead to confusion. 

The resolution of a design is the length of its shortest word. Maximum 
resolution is the usual criterion to select a design. However, for given k 
and p, and specially when these values are large, there may exist many 
different maximum resolution 2^~'p designs. For instance, the five different 
2^~^ designs are of resolution III. Additionally, for > n/2 all designs 
but one (the exception corresponds to k = n/2 and has resolution IV) are 
of resolution 

Fries and Hunter (1980) introduced the notion of aberration as a more 
powerful method to compare fractional factorial designs. For two 2^~p de- 
signs di and d2, suppose r is the smallest value such that ar{di) / ar{d2)- 
We say that di has smaller aberration than d2 if ar{di) < ar{d2)- The design 
d is of minimum aberration (MA) if there is no other design with the same 
number of factors and runs having less aberration than d. 

This previous analysis is of practical interest for designs with k > m, 
which are called fractions. The usual way of defining a fraction is the fol- 
lowing: first m initial factors b = {Bi, B2, B^} are taken to form the 2™^ 
design, and then p additional factors Bm+i, Bm+2i B^+p are assigned to 
a subset {4'i)^2) ■■■^i'p] ) selected from b = IIm\b. The selection is denoted 
by 

Bm+l = ^,l{Bl,B2, 
Bm+2 = ^2{Bl, B2, 

Bm+p = i.piBl,B2, 

The design d is thus 

d={B^,B2,...,Bm,e„^'2,...,Q (4) 

This definition requires the inclusion of the m basic columns. Any set of 

Hm containing m independent columns provides designs equivalent to those 
defined from the above criterion. If we assume that this method has been 
applied to define the 2^^"^ design, the fraction is uniquely determined for 



■,Bm) 



• ) Bm) 



(3) 
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'^2' •■•) ^p- Multiplying each relation Bm+j = by Bm+j we obtain / = 
S,[Bm+j- The p terms S^'^Bm+j are called generators of the fraction and 
represent effects confounded with the mean. These relationships are usually 
rewritten as 

I = i'lBm+l = C2-^m+2 = • • • = ^^'pBm+p- (5) 

The product of two generators of (|5|) provides another effect that is also 
confounded with the mean. The set of distinct words formed by all possible 
products involving the p generators gives the defining relation of the fraction. 
It is usual to refer to letters instead of factors or columns. A letter is any of 
the labels Bi used to denote a factor. In this setting, as before, the product 
of two identical letters is the identity. The word length is simply the number 
of letters of a word. 

For instance, if we call E, F, G, H and J the five additional factors, the 
generators of the design ([l|) are 

/ = ABCE = ABDF = ACDG = BCDH = ABCDJ. (6) 

(The non basic factors have been underlined). By taking all possible prod- 
ucts in which the five generators of (^) are present, we obtain the 32 words 
(31 if we exclude I) which form the defining relation of the fraction. From 
the lengths of the words we obtain the following WLP, 

(0,0,4,14,8,0,4,1,0). (7) 

The resolution of this design is ///, which is the maximum for 2^~^ designs. 

For the cases to be considered, k > n/2, any subset has always m inde- 
pendent columns and always leads to designs which can be defined by means 
of the scheme described in (^). In the following sections we shall use any of 
the two forms of identifying a design. 

The WLP is an essential instrument to evaluate a 2^~'p design and its 
computation requires p independent generators. In section 5 we shall en- 
counter designs which are defined for subsets of Hm that do not include the 
basic columns such as 

d2 = {A, BC, BD, CD, ABC, ABD, ACD, BCD, ABCD}. 

A simple method to obtain the generators consists in transforming d2 by 
means of an isomorphism into another set with structure (Q). For instance, 
the isomorphism T defined by T{A) = A, T{B) = BCD, T{C) = ACD and 
T[D) = ABD transforms d2 into the equivalent design 

{A, B, C, D, AB, AC, AD, BC, ABC}. 
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Using E, F, G and J as labels for the non basic factors, the generators are 

/ = ABE = ACF = ADG = BCH = ABC J. 

The WLP of d2 is (0,0,8,10,4,4,4,1,0). Since 03(^1) = 4 and 03(^2) = 8, 
di is preferred from the aberration point of view. 

Tang and Wu (1996) propose a new approach to characterize MA designs 
in terms of their complementary designs. Let Hm = dU d , they argue tliat 
when the elements in d are more "dependent", those in d should be less 
"dependent" and thus may have less aberration. They developed a general 
theory to support this intuition. Here we give a different version of Tang 
and Wu's intuition based in the concept of rank which we define as follows. 

Definition 2 The rank of a set f of columns of H^n is the number of inde- 
pendent columns it contains. 

The rank of Hm is m and hence the rank of any subset / C is never 
larger than m. Moreover, the maximum number of elements of a subset / 
with rank v is 2" — 1. Therefore, if the number of elements of a set is h, its 
minimum rank is 

Vh = \^og2{h + l)], 

where [(a)] is the smallest integer larger or equal to a. 

The maximum resolution (and also the MA) fractions always have max- 
imum rank equal to m. In the next section we shall also see that the com- 
plementary design of a MA design has minimum rank. 

Let us go back to figure 1 to identify the MA designs. For each value 
of k in the figure, we have identified with a black node the MA design. 
The complementary (symmetric) of a MA design is the worse design from 
the aberration point of view. All designs with A; > 8 are of resolution III 
with the exception of a design for k = 8 having resolution IV (the circled 
one). This design is called saturated design with resolution IV and plays an 
essential role when building minimum aberration designs. We can observe 
that all minimum aberration designs for A; > 8 start from it. This property 
and other related ones are generalized in the following sections. 

3 Screening designs {k > n/2) of minimum aberra- 
tion 

Let d be the 2*^"^ design with n/2 < k < n — 1, defined by the generators 
I = C'lBm+i = ^2^m+2 = • • • = C'pBm+p- Our goal is to choose ^'1,^2, ->Cp 
in such a way that the design d is of minimum aberration. 
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The solution for the case k = n/2 is well known. Given the basic columns 
Bi,B2, ...,Bm, it consists in assigning the (n/2 — m) non-basic factors to the 
columns obtained as the product of an odd number of basic ones. There 
is only one possible design of resolution IV and hence it is of minimum 
aberration. 

Let ipi,ip2, ■■■,(Pg with q = 2™~^ — m, be all the odd combinations of 

the basic columns Bi, B2, Bm- The generators of the 2^y^ design with 
I = 2"^-^ are 

/ = ipiBm+l = ^2Bm+2 = • • • = (f^Bm+q (8) 

From now on we will denote this design as Om- 

The extreme case is obtained when k = n — 1, and corresponds to the 
maximum number of factors which can be analyzed with a regular fraction 
of n runs. This design is called a saturated design of resolution III (or 
just a saturated design) and will be denoted as Hm, since it contains all its 
columns. 

When the number of factors to be studied is n/2 < k < n — 1, the 
maximum resolution of the design is III. The resolution criterion does not 
discriminate among designs with this number of factors and it is necessary 
to use the aberration criterion to differentiate between them. The following 
theorem gives a condition for a design with a number of factors k > n/2 to 
be of MA. 

Theorem 1 Let d C he a 2^~'p design and d its complementary, Hm = 
d\J d. A necessary condition for d to he of minimum aberration is that the 
set d has minimum rank. 

The proof, based on a combinatorial argument, is given in appendix 1. 

It should be noted that theorem 1 does not limit the number h of columns 
of d. When h > n/2, the rank of the set d of columns to be eliminated is 
always m and the theorem is trivial. The result is important for h < n/2, 
which corresponds to designs of resolution ///. We wish to emphasize the 
interest of this result when looking for MA fractions. For example, the 2^'^~^ 
design with generators 

/ = ABCE = ABDF = ACDG = BCDH = AD J = BDK = CDL = ABCDM, 

(9) 

is obtained by eliminating from the columns d = {AB, AC, BC} , which 
form a set of rank 2. The WLP of (^) is 

(0, 0, 16, 39, 48, 48, 48, 39, 16, 0, 0, 1). (10) 
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If instead of eliminating the columns indicated above we choose {AB, AC, ABCD}, 
which has rank 3, the fraction obtained has as WLP 

(0, 0, 17, 38, 44, 52, 54, 33, 12, 4, 1) 

that has more words of length 3 and is therefore worse from an aberration 
point of view. 

When the number of columns eliminated is h = 2^ — 1, where v is any 
integer, the choice of a set d with minimum rank ensures that the resulting 
fraction has minimum aberration. This is due to the fact that all sets of 
minimum rank and 2*^ — 1 elements are isomorphic, Tang and Wu (1996). 

When k > n/2, the number of generators needed to define a design is 
larger than the number of odd combinations of the basic factors. We will 
show that the minimum aberration design is obtained by adding r = k — n/2 
additional columns to Om- In this case, q out of the p generators needed are 
the same as those used in the saturated design (^). Hence, we just have to 
find the additional r = p — q generators among the 2"^"^ — 1 free columns of 
Hm, which correspond to the even products of the basic ones Bi,B2, ...,Bm- 
We shall denote these columns as ^|J^, to distinguish them from the odd ones 
(Pj-, the set they form will be denoted as 

Em = {i^i,i^2,-,ip2"^-^-i}- (11) 

Em C Hm and {I}UEm forms an abelian subgroup of {I}UHm- The sets Em 
and Hm~i are isomorphic. The problem of choosing r columns from within 
the 2™"^ — 1 even columns {Em) is the same as the initial problem (choosing 
k from Hm), though from a set of lower dimension Hm-i- Therefore, the 
problem can be solved through a iterative procedure. 

Theorem 2 Let d C Hm be a 2^~p design with n runs and k > n/2. A 
necessary condition for d to be of MA is Om C d. 

Proof. Let Hm = d U d, where the number of columns oi d is h = 
(2™ - 1) - k. Since k > n/2, then h < n/2 = 2"^"^ and, by theorem 1, 
the rank v oi d should verify v < m — 1 for d to be of MA. The set Em is 
isomorphic to any set with all its 2*"^^ — 1 elements generated by m — 1 
independent columns. Let d' be the set isomorphic to d such that d' C Em- 
The complementary of d' contains all the odd combinations of the basic 
columns and is isomorphic to d. □ 

Another equivalent statement for the above theorem in terms of the 
defining relation is the following: any minimum aberration 2^~P design d 
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with n = 2^ P runs and k > n/2 can be defined by 

/ = LpT^Bm+l = (/?2-^m.+2 = • • • = ifgBm+g = 

= Ip'lBm+q+l = 4''2Bm+q+2 = • • • = 1p'j.Bm+q+r (12) 

where ipi,ip2-, ■■■^'^q are all = n/2 — m columns obtained as odd products 
of three or more of the basic columns Bi, B2, and V'l) '02' V'r 

even products of the basic set. 

4 Saturated designs of resolution IV {Om) 

The strategy we presented in the next section to build resolution /// designs 
takes Om C Hm, the saturated designs 2^"'^ of resolution IV with generators 
(P), as the starting point. For these designs / = 2™"^^ = q + m. The smallest 
of these designs is 2^-1, followed by 2^"^, 2^^-^"" , 2^2-26 and so forth. 

The 2^-^ design defined by the columns O3 = {A, B, C, ABC} has the 
generating equation / = ABCD and WLP (0,0,0,1). The next design is 
O4, a 2^~^ fraction with generators 

/ = ABCE = ABDF = ACDG = BCDH (13) 

and WLP: (0, 0, 0, 14, 0, 0, 0, 1). 

The defining relation is one of the 2"* alias chains which appear in a 
design d C Hm- The others correspond to the 2"* — 1 columns or elements of 
Hm- Each chain includes 2'' confounded effects. The set of these confusion 
chains is called the confusion structure. The confusion structure of Om 
designs with generators given in (^), presents certain special properties of 
interest. The q generators of these fractions are even and therefore the 
2'^ words of the defining relation are even. The alias chains associated to a 
main effect are obtained by multiplying the generating equation by the single 
letter associated to the factor; the result is a chain with 2'^ odd effects. Let 
C2r+i,j be the number of effects of size 2r + 1 in the alias chains of Bj, it 
can be seen that for all r. 



C2r+l,l — C2r+1,2 — • • • — C2r+l,l — 




On the other hand, the (Z — 1) columns of Hm not included in the set Om 
are the product of an even number of basic columns ipj S Em- For the same 
reason as before, the alias chains associated to these columns are formed by 
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2'^ effects which are even. If b2r,j is the number of effects of size 2r in the 
ahas chains corresponding to the even effect ipj G Em, then for any size 2r, 



'>2r,l — 02r,2 — • • • — W2rJ-l 



(15) 



These chains, which we shah cah even chains, play a fundamental role in 
the computation of the word-length pattern for designs of resolution /// 
with generators given in (|l^). Given a2r, the value = &2r,j is obtained 
through the following argument: the total number of effects of size 2r which 
appear in the whole confusion structure of the fraction is (2^^.) , a2r of which 
are words of the defining relation and the remaining effects are distributed 
in the (l — 1) even alias chains, i.e. 

asr + a - l)62r = Qr) ' ^ = 1'2,-,V2 (16) 

Since / = 2™"^, the values a2r and 62r only depend on m and are easily 
obtained from the above expression and the following theorem. 

Theorem 3 The word length pattern, 

W = (01,02, .■.,ai), 

of Om, the saturated design 2'"'^ of resolution IV, where I = 2™^^ and 
q = 2™"^ — m, is obtained by means of the following recurrence laws: 

a2r+i = 0, r > 0; 

02 = 0, 

1 " 

0'2r+2 



2r + 2 



I . 

/ - 2r)02r 

2r + 1 , 



(r > 1). 



Proof. The design generators are even words and the product of any 
subset of them is an even word, hence in the defining relation there exists 
no even word and 02r+i = for all r. 

Since there are a2r words of size 2r, the total number of letters which 
appear in words of this size is 2ra2r- If Oi2r,j is the number of times in which 
letter Bj appears in words of size 2r, we have 

I 

y^a2r,i = 2ra2r, (17) 
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and according to (14) and (|T^), a2r,i = a2r,2 



= Oi2r,i, for all r, therefore 



2r 

a2r,j — —a2r- 

and the number of words of the same size 2r which do not contain the letter 

Bj is a2r — Oi2r,j = (1 " 2r/l)a2r- 

To obtain the confusion chain corresponding to Bj, we multiply each 
word of the generating relation by that factor. All the effects obtained will 
be odd. The effects of size 2r + 1 in this chain correspond to either words 
of size 2r in the defining relation which do not include factor Bj or words 
of size 2r + 2 in the defining relation which do include factor Bj. If C2r+i,j 
is the number of effects of size 2r + 1 in the chain, we have 

C2r+l,j = {a-2r — "2r,j) + OL2r+2,j 

l-2r 2r + 2 

— ^ a2r H ^ a2r+2- 

All these odd effects can be found in the confusion chains of the main effects. 
The total number of effects of size 2r + 1 is Y\\-^ C2r+i 7 = I „ 1 and 
therefore 

[l - 2r) a2r + (2r + 2)a2r+2 = 

Solving for a2r+2 we obtain the recurrence law. Since the design is of reso- 
lution IV, the initial value of the law is 02 = O.D 

5 Word-length pattern 

We will show a simple way of obtaining the WLP of a 2^""^ design with 
generators (p!2[). For this purpose we shall now introduce the concept of 
word-length pattern polynomial (WLPP). 

Definition 3 Let d C Hm be a 2^^^^ design with WLP, W = (oi, 02, Ofe). 
The polynomial P^ associated to this design is given by 

Pd{u) = 1 + aiu + a2U^ + • • • + a^u^ . 

We shall refer to it as the word length pattern polynomial (WLPP). 

The above definition is valid for the case of full designs where Pd{u) = 1- 



I 

2r + 1 
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For simplicity we shall denote by = Pq^ the polynomial associated 
to Om, whose coefficients a2r can be obtained by applying Theorem |3|. Let 
us call 

2m-l 

Qm{u) = ^ b2rU'' (18) 
r=l 

the polynomial with coefficients b2r = b2r,j defined in the preceding section. 
The equation (16) can be rewritten as 

^ /om-l\ 

Pm{u) + {2^-'-l)Qm{n)=Y,( 2r r' ^^^^ 

We will refer to Qm as the effect-length pattern polynomial for the even alias 
chain. 

Let d C Hm be the 2^~p design defined as d = Om U e, where Om 
is the design (P) and e C Em with r = k — n/2 columns. Consider the 
design e, a design formed with columns from Em- The set Em is isomorphic 
to Hm^i, therefore any design e has a maximum of m — 1 independent 
columns and is isomorphic to a design of Hm-i- The word length pattern of 
e can be obtained taking into account the definition of word given in (^) or, 
alternatively, it could be find the isomorphic set of Hm-i and write it as in 
(§). For example, the design with generators (P) contains both O4 and the 
subset 

e = {AD, BD, CD, ABCD} 

of 

Ei = {AB, AC, AD, BC, BD, CD, ABCD} . 

E4^ is formed by all the possible even combinations of the basic columns 
A, B,C, D. It can be seen that the rank of E4 is 3. Using the additional 
factors in ^ to define the generators J = AD, K = BD, L = CD and 
M = ABCD, the set e is isomorphic to {J,K,L,JKL}, which defines a 
2'*"^ design of H^. The defining relation of e is / = JKLM and hence its 
WLPP is Pe{u) = 1 + 7x1 

Theorem 4 Let d be a 2'^"^ design defined by d = Om U e, where Om C Hm 
is the saturated design of resolution IV and e C Em- The (WLPP) of d is 
given by 

Pd{u) = (1 + uYQm{u) + Pe{u)[Pm{u) - Qm{u)], (20) 
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where r = k — n/2, Pm and Pe are the WLPP of Om and e, respectively, and 
Qm{u) = 

Proof. Suppose that the design d has only one more column than Om- 
In this case, the generators of d are the q = n/2 — m generators of Om given 
in (P), plus Bm+q+iipi, where V'l has an even number of basic factors. The 
defining relation has 2"?+^ words, the first 2"? correspond to the design Om 
which has WLP polynomial Pm- The remaining 2'^ words are the result of 
multiplying the 2'? words of the defining relation of Om by Bm+q+itp'i- When 
multiplying the defining relation of Om by Bm+q+i^Ji we obtain the alias 
chain associated to ip'i plus the new letter Bm+q+i in each effect. The effects 
of size 2i in the alias chain of ipi become words of size 2i + 1 in the defining 
relation of d. Since the effect length pattern polynomial of the alias chain of 
tp'i is Qm{u), it will be uQm{u) for the chain Bm+q+iipi so that 

Pd{u) = Pm{u) + uQm{u)- 

Assume now that d is given by Om plus two additional columns, i.e., we 
add Bm+q+itpi and Bm+q+2ip2 to the generators of Om- The defining relation 
of d is formed by the defining relation of Om plus the words resulting from 
multiplying the defining relation of Om by Bm+q+iipi, Bm+q+2i^2 and those 
corresponding to the product of both, Bm+q+iBm+q+2'>Pi''p2- The product 
ipiip2 is an element of Em different from V'l and ^^2- The contributions of 
each new generator to the generating equation of d will be: uQmiu), corre- 
sponding to Bm+q+iip'i', the same value uQm{u), corresponding to Bm+q+2ip'2 
and v^Qmiu) for Bm+q+iBm+q+2'ipi'ip2- The term is due to the fact that 
all the effects of the alias chain V''iV'2 ^i^^ contain the two new letters 
Bm+q+iBm+q+2- The WLP polynomial of d is 

Pd{u) = Pm{u) + 2uQm{u) + V?Qm{u) 

adding and substracting Qm{u) to the right hand side, we obtain 

Pd{u) = (1 + ufQmiu) + [Pm{u) - Qmiu)]- 

The above result is generalized when d = Om U e, and the elements of 
e = {V'l, '02) •■•> V'r} are independent {r < m — 1). The independence of 
the elements of e ensures that any product of a subset of them provides a 
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different element of Em. With r generators one can get combinations of 
i of them that will incorporate {\)u'^Qm{u) to the polynomial Pd{u), 



Pd{u) = (21) 
i=l 

= {l + uYQm{u) + [Pm{u)-Qm{u)] 

When the elements of e = j^'j^, ^'2, V'r} ^^^^ independent (this is 
the case where r > m — 1), the product of some subset of them is equal to 
/ instead of some element of Em- Thus, if 

the product of the generators Bm+q+nip'j^, Bm+q+jii^'j^, Bm+q+j.ip'j^ is 

Bm+q+ji Bm+q+j2 ' ' ' ^m+q+ji 

and its contribution to Pd{u) is u^Pmiu) instead of u^Qmiu)- Let be the 
number of subsets of i elements of e the product of which is equal to /. The 
polynomial Pd{u) is 

Pd{u) = PM + ^[.ju'Qmiu) + Y,ay[Pm{n) - Qm{u)], (22) 

where the two first terms from (21) correspond to the counting process if the 
r elements of e are independent, and the term Yll=i'^i^^[Prn{u) — Qmiu)] 
corrects this value by taking into account the dependencies on e. Once again, 
adding and substracting Qm (u) to the right-hand side of ( [2^ ) and regrouping 
terms one gets 

r 

Pd{u) = (1 + uYQm{u) + (1 + a[u')[Pm{u) - Qra{u)]. 



1=1 



The polynomial Pe{u) = (1 + X^[=i a^^*) is by definition the WLPP of design 
e. This last result is the expression ( pO[ ) which we wanted to prove. □ 

Example 2. The WLP for the 2^'^"^ minimum aberration fraction with 
generators (0) is obtained by applying (E 



Pd{u) = (1 + uYQi{u) + P,{u) [P^{u) - Qi{u)] 

where Pi{u) = 1 + 14n^ + is the polynomial corresponding to the fraction 
(ll), and Qi{u) = Au"^ + + 4u^ is obtained from (|l|) and (|l|). Also, 
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Pe{u) = 1 + is the polynomial of the 2^ ^ minimum aberration design 
(/ = JKLM). By substitution in we obtain: 

Prf(u) = 1 + iQu^ + 39n^ + 48u^ + 48ti^ + 48n^ + 39n^ + 16ti^ + 

The coefficients of Pd are the elements of the word length pattern (p!o|). 

Using theorem it is easy to prove that the goodness of the design d in 
terms of its aberration depends entirely on e. 

Theorem 5 Let d he a design defined by d = Om U e, where Om C is 
the saturated design of resolution IV and e C Em- The design d is MA if 
and only if e is isomorphic to an MA design e' C Hm-i- 

Proof. Let di C Hm and d2 C Hm be two different 2'^"^ fractions. The 
design di is of smaller aberration than d2 if and only if (by definition of 
WLPP) Pdi{u) — -Prfal^) a negative coefficient for the term of lowest 
order. If di = Om U ei and d2 = Om U 62, then by (pO|) we have 

Pd,{u) - Pd,{u) = [PeAy) - Pe,{uWm{u) " Qm{u)] 

and since Pm{u) — Qm{u) = 1 — 62"^^ + - • • , if the coefficient of the lowest order 
term of Pei {u) — Pe2 (u) is negative, the corresponding one for Pd^ (u) — Pd2 (u) 
will also be negative. Therefore, if ei is isomorphic to an MA design of Hm-i, 
di will also be MA, and vice versa. □ 

This theorem is useful to justify that the design d with generators (^) is 
MA. This follows from the fact that this design is the union of O4 and e = 
{AD, BD, CD, ABCD} , and the subset e is isomorphic to {J, K, L, JKL}, 
which represents the 2'^~^ MA design. 

Knowing the 2^~^ MA (m = 3) design it is possible to obtain the 2^^~® 
MA (m = 4) design. This design in turn allows the determination of the 
228-23 y^j^ _ 5^ design, and so forth. Given the MA design d C Hm 
with k factors, one can immediately obtain the MA design d' C Hm+i with 
2™ + A; factors. 

The d C Hm MA design with a number of factors k > n/2, is obtained 
as d = Om U e, where Om is the (|8|) design and e C Em with r = k — 
n/2 columns. To determine e one should obtain the MA with r factors in 
Hm-i- Let d' C Hm-i be the MA design with r factors defined from the 
basic factors B[, B2, B!^_i. The set e is defined as the transform of d' 
by the isomorphism T from Hm-i to Em defined by T{B'j) = BjBm for 
j = 1,2, ...,m — 1. When r > n/4 the method can be applied iteratively. 

Example 3. Let us illustrate, for instance, how to obtain the 2^^~^^ 
MA design. This design is included in i?5, and, following theorem 2, its 
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first 16 columns correspond to all the odd combinations (O5) of the basic 
columns A, B,C, D, E. The remaining 12 columns arc a subset of the even 
combinations of the basic ones (£^5) isomorphic to the 2^^^^ MA design 
included in if4. The 2^^~^ MA design (verifying 12 > n/2 = 8) should 
contain the 8 odd combinations of the 4 independent columns plus 4 other 
combinations of them. We take cq = {AE, BE, CE, DE} as the independent 
columns of E5 and build with them the 8 odd combinations, 

ei = {AE, BE, CE, DE, ABCE, ABDE, ACDE, BCDE}. 

The 7 even combinations of cq form the set {AB, AC, AD, BC, BD, CD, ABCD} 
isomorphic to H3. Prom it we take 

62 = {AD, BD, CD, ABCD} 

which is isomorphic to the 2^~^ MA design of H^. The 2^^~^^ MA design is 
d = O5 Uei Ue2. 

The case 2^^~^^ has been chosen to illustrate the proposed method. A 
simpler way to build the design in this case is to eliminate 3 columns out 
of the 31 which form H^. By results of section 2, it suffices to choose 3 
minimum rank columns of H^,. In this case, d is missing just the columns 
d = {AB,AC,BC} to complete H^, and the rank of d is the minimum 
possible. 

6 Conclusions 

The selection of a 2''-P re gular fraction implies choosing k columns among 
the n = 2"* available in a full factorial design with m = k — p factors. When 
k > n/2, it can be shown that the minimum aberration design contains n/2 
columns obtained as a product of an even subset of m independent columns. 
The choice of the other k — n/2 columns requires solving the same problem 
but now in a smaller set, the full factorial design with m — 1 factors. A 
iterative procedure has been derived from this property that allows building 
screening designs {k > n/2) of minimum aberration for large k. 

The comparison of designs in terms of aberration requires the computa- 
tion of word-length patterns. When k and p are large the direct counting 
process is very cumbersome. In this article we introduce a simple relation- 
ship, based on products of polynomials, for computing the word-length pat- 
tern for these designs. This relation provides an interesting correspondence 
between designs of different number of runs. 
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7 Appendix: proof of theorem 1. 



We choose d = ^2) ChJ so that the complementary design d has the 
minimum number of words of size 3 in its defining relation. In the complete 
defining relation of Hm, any factor (letter) appears in (n — 2)/2 words of 
size 3, therefore, if we remove h columns, the number of words of size 3 that 
we eliminate is 

_ n — 2 — — 

a{d) = {^-)h-a'^{d)-2a3{d), (23) 

where a'^{d) and 03 (d) are, respectively, the number of these words which 
contain two and three eliminated factors. Moreover, in the complete gener- 
ating equation of Hm any combinations of the two factors appear once in 
words of size 3, thus 

(2) =«3®+3a3(d). 

Substituting in (^) we find that the number of words of size 3 that have 
been eliminated is 

h(n - h-1) .-. 
a{d) = + 03(0!). 

and this number will be maximum whenever 03 (d) is also as large as possible. 
The term a^{d) denotes the number of relations or words 

UX = I, ^-^thl,~Ct<^d. (24) 

We shall now see that if 03 (d) is maximum, then d is generated by the 
smallest number of independent columns Vh S N, with h < 2'"'^ (note that it 
is not possible to generate h different columns using less than independent 
columns) . 

Suppose that the number of independent elements of (i is w > w/i, we 
shall now see that it is then possible to define a new set with h columns 
which provides a larger number of words of size 3, which contradicts the 
statement that a2,{d) is maximum. 

Let Jy be the set formed by 2^^ — 1 different columns generated by the v 
independent columns of d. Out of the v independent columns we choose f — 1, 
which we will call tpi,ip2, tpy-i- Let Jv~i C Jy be the subset generated by 
these columns and Jv-i the complementary set, such that Jv-i = Jv\Jv-i- 
We divide the set do = d into two parts, 

eo = do n Jy-i = {^1, •••> V's} 
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and 

For this proof we have used on purpose the notation tp^, ipj for the columns 
of Jy, Jv-i, respectively, which, without loss of generality, can be considered 
as even and odd combinations of the basic factors. As will be seen below, 
this notation helps to understand the following development, if we bear in 
mind properties such as that the product of two even combinations is also 
even and so forth. 

The elements of /o C Jv are obtained as a product of the generators of 
Jv-i V'l) ■02) ■■■,tpv-i any other term of /o, for example ipi, thus 

iPf- = (piWt-, con tiJt G Jt,_i, para i = 2, /i — s. (25) 

With the new notation, the relations (|2^ ) can be divided into two classes: 

(A) : tpiipjijf, = I, 'i/'i,V'i,V'fc e eo (26) 

and 

(B) ■■ fi'Pj'iPk = I fi,'Pj^fo, V'fceeo. (27) 

(Note that taking into account the even/odd features no other combinations 
are possible). 

We now replace € /o by ^p^ G Jv-i (V't ^ ^o) as follows: we choose 
two different ipa^i^b ^ such that V'l = V'aV'fe ^ ^o, (since s < 2^^^ this 
choice is always possible). In (p5|), is replaced by Tp'^ and we obtain the 
columns V't S Jv-i, such that 

'4>[ = 7p[zut ,t = 2,3, ...,h - s. 

We define Cq = : ip^ = ip'iWt ^ eo}, and the new sets ei = eo Ue'o U 
and 

/i = Wt = fi'^t G /o : 01 = i^'iWt G eo, t = 2, 3, h - s}. 

If di = ei U /i, we now show that 03(^1) > 03(^0)- Since eo C ei, the other 
words of type ( p6| ) initially formed can still be formed after the substitution 
process. For the words of type (p7|), there are three possible situations: 

1. (p^, ifj G di, that is ip^, have not no been replaced in which case the 
word Lp^ifjipj. = / is maintained. 
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2. Both ^Pi,(pj ^ di, have been replaced by ip'^ = ipiWi and ip'j = ipiZUj, 
respectively. In this case ip'^ip'jipi^ = I will be a new word . Given that 
ifiipjipf. = I, according to (p5|) WiWjipj^ = /, therefore {ipi'coi) {ipiWj) ipj^ = 

3. ^ di, ipj € di,ipi have been replaced by tp[ = tpiZUi, but ipj has not, 
because the substitute V'j already existed, i.e. Tp'j = ipiWj G d^. Using 
the same reasoning as above, since (/9j(/9j'i/'fc = ^i tben V'iV'jV'fc = ^ "^i^^ 
be a word that did not exist before (given that ip^ is a new element). 

In this case it is possible that two different words give equal transforms. 
Suppose that ip^,ip^,tpg G fo:'4'ui'^v ^ ^Oj that 

^i^ri^u = I and fifsi'v = I: (28) 

and ^ di has been replaced by ip[ S ei, while ^Pr,Vs ^ /i could 
not be replaced. Moreover, it holds that ip'^ = ipiWu = ip^ and ip'^ = 
ipiWv = tpu- In this situation, the two last words ( p8[ ) have the same 
transforms ip[ip^ip^ = I. But in this case the new word ip^ip^ip'^^ = I is 
created, which did not exist since ip^ ^ do. 

A third coincidence could appear, if z = 1 and V'a = i^u and jpj, = ip^. 
In this case the transform of (^) coincides with ip'^ip^ip^^ = /, the 
relation used to define ip'i. 

The substitution process verifies that 03(^0) < 03(^1). If there is a triple 
coincidence, as described in point 3 of the analysis above, then 03(^0) < 
03(^1), if not, a3((io) < 03(^1) (all words have been transformed and addi- 
tionally Ip'i'lpa'^b = I)- 

Repeating the process, a finite sequence of sets do,di,d2, ---jdM {M < 
h — q), is created, which verify 

asido) < a3('^i) < • • • < as(dM~i) < a^iidu) (29) 

where dM C Jv-i {fM = 0)- Moreover, in the last step, a^{dM~i) < 03('^A/) 
since the situation described in point 3 can not arise, since Jm = (all 
terms of type (p^ have been replaced ). 

The final set dM obtained has more words than d, which contradicts the 
statement that as{d) is maximum. Therefore, the theorem is proved. □ 
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